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ABSTRACT 


It  is  wall  established  that  tha  rata  of  tharaal  oxidation  of 
silicon  depands  on  the  surface  orientation,  but  a  comprehensive  model 
for  the  role  of  the  surface  orientation  in  the  kinetic  mechanism  is 
lacking.  The  results  of  several  experiments  designed  to  develop  a 
better  understanding  of  which  surface  properties  are  most  important  in 
establishing  the  oxidation  rate  are  reported.  The  results  indicate 
that  the  Si  surface  atom  density  is  important  in  the  initial  stages  of 
oxidation  and  a  revised  explanation  for  the  crossover  effect  is 
proposed . 


Introduction 


The  silicon  surface  orientation  is  an  important  parameter  in  the 
kinetic  mechanism  for  the  thermal  oxidation  of  Si  (1-11).  Orientation 
dependent  oxidation  rates  are  observed  at  all  temperatures  in  the 
range  SOO-llOO^C  (1-11).  The  order  of  the  rates  as  a  function  of 
orientation  is  complex  and  depends  on  oxide  thickness  (7-9) ,  oxidation 
temperature  (3,7-9),  and  ambient  pressure  (1,6-8).  In  the  initial 
stages  of  oxidation  at  1  atm  the  rates  are  in  the  order 

(110)  >  (111)  >  (100),  but  within  the  first  25  nm  of  oxide  formed,  the 
order  is  reversed  on  the  (110)  and  (111)  orientations  at  temperatures 
in  the  range  750- 1000,:aC  (7,9-11).  The  thickness  where  the  rate 
crossover  occurs  varies  slightly  with  oxidation  temperature, 
increasing  with  increasing  temperature  (7) .  Further  complicating  an 
understanding  of  the  role  of  orientation  in  the  kinetic  mechanism  is 
the  crossover  observed  between  the  (111)  and  (lOO)  orientations  at 
partial  pressures  below  0.1  atm  (6-8).  None  of  the  current  models  for 
Si  oxidation  can  explain  all  aspects  of  the  observed  orientation 
dependent  rate  behavior  described  above.  In  order  to  model  the  effect 
of  surface  orientation  on  the  oxidation  kinetics  it  is  necessary  to 
identify  which  surface  properties  alter  the  rate  and  which  steps  in 
the  kinetic  mechanism  are  affected. 

The  oxidation  process  involves  transport  of  the  oxidant  species 
to  the  Si-Si0e  interface  followed  by  a  chemical  reaction  at  the 
interface  (12,13).  The  most  widely  accepted  model  for  Si  oxidation  is 
the  Deal-Grove  model  (14)  which  uses  a  linear-parabolic 
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expression  (15)  to  fit  the  thickness,  L,  versus  time,  t,  data  aa 
followa: 

<L*  -  Lo*>/  kp  ♦  <L  -  Le> /  kx  *  t  -  t«. 

Thia  aodal  ia  darivad  froa  a  ataady  atata  analyaia  of  a  diffuaive 
tranaport  flux  in  aariaa  with  a  firat  ordar  intarfaca  reaction.  The 
parabolic  rate  conatant,  kp,  ia  proportional  to  the  effective 
diffuaivity  of  oxidant  through  the  SiO«  film,  and  the  equilibrium 
concentration  of  oxidant  in  the  oxide,  C-.  The  linear  rate  conatant, 
ki,  ia  proportional  to  the  interface  reaction  rate  conatant,  k,-,  and 
C~.  L0  and  t0  are  offsets  in  thicknaaa  and  tine,  respectively,  and 
define  an  initial  regiae  where  the  oxidation  rate  ia  greater  than 
predicted  by  the  model.  Applying  thia  model  to  experimental  data  it  ia 
found  that  both  kp  and  kv  are  orientation  dependent  (2, 4, 5, 9). 

Although  the  Deal -Grove  model  provides  a  basis  for  analyzing  the 
kinetic  data  and  allows  for  orientation  dependent  rates  through  k^,  it 
does  not  explicitly  account  for  the  observed  orientation  behavior. 

Other  studies  have  suggested  possible  sources  for  the  orientation 
dependent  rates  including  variations  in  the  number  of  Si-Si  bonds 
available  for  reaction  <1,5>,  the  orientation  of  bonds  on  the 
surfaces  (1,16),  and  the  density  of  surface  steps  (17,18).  At 
temperatures  below  1000~C  the  oxidation  rate  may  be  further  influenced 
by  intrinsic  oxide  stress  (10,19) .  The  presence  of  an  intrinsic  stress 
having  a  magnitude  on  the  order  of  10*  dynes/cm“  at  oxidation 
temperatures  below  approximately  1000°C  is  well  established  (20-24). 

(2) 


Several  Modifications  to  tha  Daal-Grova  Modal  hava  introducad 
oriantation  dapandant  paraaatara  such  aa  strata  (25-28)  and  tha  Si 
surface  atom  density  (25)  into  tha  rate  constants  to  explain  tha 
orientation  dependence.  However,  it  is  not  clear  whether  stress 
increases  the  oxidation  rate  by  enhancing  the  interface 
reaction  rate  (25),  or  decreases  the  oxidation  rate  by  reducing  the 
diffusivity  of  oxidant  through  the  oxide  (26-28) .  The  goal  of  this 
study  was  to  deteraine  the  effect  of  stress  and  Si  surface  atoM 
density  on  the  oxidation  rate  at  1  atm.  There  is  not  sufficient  data 
available  on  the  variation  of  surface  steps  with  orientation,  or  on 
the  effect  of  bond  orientation  on  the  oxidation  rate  to  determine 
their  role  in  the  orientation  dependence. 

The  results  of  oxidation  experiments  at  700oC,  1000°C  and  HOO-^C 
using  five  Si  orientations  are  presented.  The  use  of  an  extended  range 
of  orientations  provided  a  broad  basis  for  comparison  of  the  oxidation 
rates  with  Si  surface  properties.  It  is  difficult  to  extract  accurate 
rate  constants  from  the  kinetic  data,  especially  at  low  temperatures 
where  the  oxide  thickness  range  covered  in  a  typical  oxidation 
experiment  is  small.  Therefore,  in  this  study  the  kinetic  data  was 
analyzed  in  terms  of  oxidation  rate  as  a  function  of  oxide  thickness, 
i.e.  the  slope  of  the  thickness  versus  time  curve  at  a  given  film 
thickness.  This  enabled  a  model  independent  analysis  of  the  oxidation 
process  as  a  function  of  both  oxide  thickness  and  Si  orientation.  For 
oxidations  at  700°C  it  was  anticipated  that  the  slow  growth  rate  would 
limit  the  oxidation  to  a  thickness  regime  where  only  the  interface 
reaction  kinetics  were  important  in  controlling  the  oxidation  rate.  At 


this  temperature  the  intrinsic  oxide  stress  ia  high  (23,24)  and  tha 
affact  of  oxida  atraaa  on  tha  intarfaca  raaction  waa  axpactad  to  ba 
aaxiaizad.  In  ordar  to  aaparata  atraaa  af facts  from  othar  oriantation 
dapandant  paraaatara,  a  aacond  aat  of  oxidation  axpariaanta  using  all 
five  orientations  waa  carried  out  at  1000°C.  Stress  relaxation  is 
observed  at  teaperaturea  above  970°C  (20)  and  thus,  at  1000SC  the 
intrinsic  stress  is  greatly  reduced  (21,23,24).  However,  experimental 
streaa  measurements  show  a  low  residual  intrinsic  stress  is  present 
even  in  oxides  grown  at  1000°C  (24).  To  determine  if  the  crossover  was 
affected  by  stress,  oxidations  were  done  at  1100°C  which  is  well  above 
the  temperature  required  for  viscous  relaxation.  The  results  lend  new 
insight  into  the  crossover  effect  and  the  importance  of  surface  atom 
density  in  establishing  the  oxidation  rate. 

Experimental 

Commercially  available  single  crystal  Si  wafers  with  a 
resistivity  of  2-10  ohm-cm  and  having  the  orientations  (110),  (111), 

(311),  (511)  and  (100)  were  oxidized  at  1  atm  at  TOO^C,  lOOO^C  and 

llOO^C  in  pure  dry  oxygen.  The  (311)  and  (511)  orientations  were  N- 
type  wafers  while  the  other  orientations  were  P-type  but  in  this 
resistivity  range  the  oxidation  rate  is  independent  of  dopant 
type  (5) .  All  oxide  thickness  measurements  were  made  with  a  research 
quality  ellipsometer  using  632.8  nm  light.  The  samples  were  cleaned 
immediately  before  oxidation  using  a  modified  RCA  (29)  procedure 
followed  by  a  10  sec  HF  dip  and  DI  water  rinse.  It  has  been 

(O 


established  that  the  initial  oxidation  rataa  vary  with  tha  wafar 
cleaning  procedure  (30,31).  To  inaure  the  validity  of  the  reaulta 
preaented  here,  all  aaaplea  were  cleaned  identically  and  the 
concluaions  are  baaed  on  relative  rather  than  abaolute  oxidation 
ratea . 

The  700,=C  oxidation  kinetica  were  meaaured  in-aitu  with  an 
oxidation  furnace  mounted  at  the  focua  of  the  ellipaometer  (32). 
Optical  conatanta  were  obtained  from  previoualy  publiahed 
data  (22,32).  The  samplea  were  heated  to  the  oxidation  temperature  in 
N^;  for  15  min  before  changing  to  an  0a  ambient.  The  0**  and  N*.  were  of 
high  purity  having  an  Ha0  content  of  leaa  than  1  ppm  meaaured  at  the 
outlet  of  the  oxidation  furnace.  Kinetic  data  waa  collected  over  46-52 
hour  time  intervala.  The  oxide  thickneaa  obtained  in  thia  time  waa 
less  than  30  nm  on  all  samplea.  The  oxidation  rate  was  alow  enough 
that  manual  null  el 1 ipaometr ic  measurements  could  be  made  during  the 
oxidation  to  obtain  accurate  thickness  values  (approximately  IX 
error )  . 

The  kinetic  data  at  lOOO^C  waa  not  collected  in-situ  but  instead 
from  a  series  of  oxidations  using  a  conventional  oxidation  furnace 
having  a  double-wa 1 led  fused  ailica  furnace  tube.  A  preburner- 
condenser  apparatus  (5)  was  placed  in  aeries  with  the  oxidation 
furnace  to  remove  H^O  and  hydrocarbons  from  the  oxidizing  ambient.  The 
measured  H^O  content  was  less  than  3  ppm  at  the  furnace  outlet. 

Kinetic  data  was  collected  from  eight  oxidations  ranging  in  time  from 
0.5  to  16  houra  and  covering  the  thickness  range  25-300  nm.  Samplea  of 


the  five  orientations  wara  oxidized  aiaultanaoualy ,  raaoved  from  tha 
furnaca  and  aaaaurad  allipaoaatrically .  For  aach  oxidation  run  the 
aaaplaa  wara  flushad  with  dry  N*  in  tha  and  cap  then  brought  to 
oxidation  temperature  at  tha  canter  of  the  furnace  hot  zone  in  an  N* 
aabiant.  Tha  aaaplaa  alao  received  a  15  minute  post  oxidation  anneal 
in  Na.  The  thickneaa  data  ia  the  average  of  four  el 1 ipsometr ic 
meaauraaenta  (2  per  wafer  with  2  wafera  of  each  orientation).  This 
type  of  experiment  confirma  the  order  of  the  relative  rates  since  all 
orientationa  are  oxidized  aimul taneoualy .  In  the  in-situ  experiments 
each  orientation  ia  oxidized  in  aeparate  runs,  end  thus,  the  relative 
rates  may  vary  slightly  aince  it  is  difficult  to  maintain  the  exact 
oxidation  conditions  from  run  to  run. 

At  HOO^C  extensive  data  aeta  were  collected  on  the  (111)  and 
(110)  orientationa  and  two  oxidatione  included  all  five  orientations. 
The  aame  apparatus  described  for  the  lOOO^C  oxidations  was  used  at 
HOO^C,  but  the  pre-oxidation  heating  and  post-oxidation  annealing 
were  done  in  an  Ar  ambient  to  prevent  the  formation  of  an  interfering 
SiaN*  film  (33).  The  H«0  content  in  the  oxidizing  ambient  was  less 
than  4  ppm  measured  at  the  furnace  outlet. 

Results 

As  explained  above,  the  data  was  analyzed  in  terms  of  oxidation 
rate  as  a  function  of  thickness  and  orientation.  The  results  of  the 
700^0  oxidation  experiments  plotted  as  oxide  thickness  versus  time  are 


shown  in  Figure  1.  There  is  a  distinct  orientation  dependent  oxidation 
rate  on  all  five  orientations.  The  initial  oxidation  rate  decreased 
rapidly  with  thickness  on  all  orientations  and  could  be  fit  to  two 
decaying  exponential  functions,  with  one  extending  to  about  3  ns  and 
the  second  extending  from  3  to  10  nm.  Similar  rate  behavior  in  the 
initial  stages  of  oxidation  was  reported  by  Hassoud  et .  al.  for 
oxidations  in  the  temperature  range  800-1000°C  (7,8).  Beyond  10  nm  the 
oxidation  rate  as  a  function  of  thickness  was  determined  analytically 
from  a  second  degree  polynomial  fit  to  the  data.  This  corresponds  to  a 
linear-parabolic  form  for  the  equation.  Oxidation  rates,  i.e.  the 
slope  of  the  fitted  functions,  at  700°C  for  several  oxide  thicknesses 
on  the  five  orientations  are  presented  in  Table  1.  The  oxidation  rates 
are  in  the  order: 


(110)  >  (111)  >  (311)  >  (511)  >  (100) 

up  to  a  thickness  of  25  nm .  However,  the  difference  in  rates  on  the 
(110)  and  (111)  surfaces  decreases  with  thickness  and  a  crossover  in 
rate  to  (111)  >  (110)  occurs  near  25  nm .  The  thickness  where  the 
crossover  in  rate  is  observed  is  greater  than  expected  based  on  the 
results  of  earlier  studies  over  the  temperature  range  750- 
1000’C  (7,11).  The  discrepancy  may  reflect  subtle  differences  in  the 
exact  pre-oxidation  cleaning  procedure  (30,31)  or  oxidation 
conditions.  A  second  crossover  thickness  has  been  defined  at  the 
actual  thickness  where  a  crossover  is  observed  in  plots  of  thickness 
versus  time  (7).  This  crossover  point  is  not  observed  in  the  700°C 


data  presented  here,  (Figure  1),  because  the  oxidation  rate  was  so 


•low  that  th«  thicknaaa  whara  a  crossover  in  L-t  occurs  was  not 


raachad  within  tha  48  hour  oxidation.  The  other  three  orientationa 
maintain  tha  saia  order  throughout  the  thickness  range  covered  in  this 
time  interval . 

The  lOOO^C  data  showed  the  same  orientation  dependence  as  the 
700,=,C  oxidations.  The  (110)  and  (111)  surfaces  had  the  highest 
oxidation  rates  with  the  initial  order  (110)  >  (111).  Because  of  the 
faster  oxidation  rate  at  lOOO^C,  the  thickness  where  a  crossover  in 
both  rate  versus  thickness  (L=10  nm)  and  thickness  versus  time  (L=40 
nm)  occurs  was  reached  within  the  first  hour  of  oxidation.  These 
crossover  thicknesses  are  consistent  with  previously  reported 
data  (7).  The  rates  on  the  (311),  (511)  and  (100)  surfaces  were  in  the 
same  order  as  at  700-=>C  with  the  thickness  range  extended  to  250  nm  at 
1000°C.  Thus,  the  order  of  the  rates  on  these  three  orientations  was 
not  changed  by  decreasing  the  intrinsic  oxide  stress.  The  data 
collected  in  the  region  before  the  rate  crossover  is  not  sufficient  to 
determine  the  rates  from  functional  fits,  but  second  degree  polynomial 
fits  were  used  to  determine  the  rates  at  thicknesses  beyond  25  nm.  The 
rates  at  lOOO'^C  relative  to  the  (100)  orientation  in  several  thickness 
ranges  are  included  in  Table  5  which  is  discussed  later.  The  order  of 
the  rates  beyond  the  crossover,  determined  in  this  study  for  dry  0e 
oxidations,  is  the  same  as  the  order  reported  by  Pliskin  (2)  for  1  atm 
steam  oxidations  at  971,=,C  using  the  (111),  (110),  (311)  and  (100) 

orientations.  However,  both  1  atm  studies  exhibit  an  orientation 
dependence  which  differs  from  the  high  pressure  results  obtained  by 


Ligenza  ( 1 ) . 


At  1100°C  the  cronovtr  behavior  on  tha  (110)  and  (ill)  aurfacaa 


waa  diffarant  than  at  lower  taaparaturaa  (700-1000*0 .  A  partial  data 
set  from  tha  oxidations  at  1100*C  on  tha  (110)  and  (111)  aurfacaa  la 
ahown  In  Table  2.  The  reported  thickneaaea  are  tha  average  of  four 
ellipaometric  measurements  on  two  samples  of  each  orientation.  The 
thickness  for  the  one  hour  oxidation  was  sufficient  to  calculate  an 
index  of  refraction  for  these  films  of  1.459.  It  is  seen  from  Table  2 
that  for  oxidation  times  extending  to  24  hours  the  oxide  thickness  on 
the  (110)  surface  is  always  greater  than  on  the  (111)  surface.  The 
maximum  variation  in  the  measured  thicknesses  was  ±1  nm  and  thus,  the 
observed  order  is  not  due  to  differences  in  thickness  across  the  wafer 
or  differences  In  the  position  of  the  wafers  in  the  furnace.  A 
crossover  is  not  observed  in  a  plot  of  thickness  versus  time  for 
thicknesses  extending  to  700  nm .  In  contrast,  for  oxidations  in  the 
temperature  range  800-1000°C  the  crossover  in  thickness  versus  time 
plots  to  (111)  >  (110)  occurs  within  the  thickness  range  15-40  nm  (8). 
The  oxidation  rates,  determined  from  second  degree  polynomial  fits  to 
the  entire  data  set,  converge  with  increasing  oxide  thickness  but  do 
not  exhibit  a  crossover.  Oxidations  on  all  five  orientations  showed 
that  the  inial  order  of  the  rates  observed  at  700=C  and  lOOO^C  is 
maintained  to  a  much  larger  thickness  (several  hundred  nm)  at  1100,=>C. 

In  summary,  the  initial  oxidation  rates  follow  the  order 
(110)  >  (111)  >  (311)  >  (511)  >  (100).  This  order  changes  within  the 
first  30  nm  of  oxide  formed  for  700,=C  <  T.=  „  <  1000,=>C  to  (111)  >  (110) 
with  the  other  orientations  retaining  the  same  order.  At  1100,=:,C  the 
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order  for  all  thicknasaaa  0-700  na  ia  Identical  to  that  obaerved  at 
lower  teaperaturea  in  the  precroaaover  region. 


Diacuaaion 

In  order  to  understand  the  role  of  the  Si  surface  orientation  in 
the  oxidation  aechanism,  it  is  iaportant  to  establish  whether  the  rate 
is  controlled  by  diffusional  transport  or  by  the  interface  reaction. 
The  rate  liaiting  process  in  these  experiaents  was  determined  by 
comparing  the  flux  of  0*  that  reacted  at  the  interface,  F„,  with  the 
calculated  diffusive  flux  of  0^,  FD,  at  the  oxidation 

temperature  (34).  The  reaction  flux  is  observed  experimentally  and  is 
calculated  from  the  oxidation  rate  as  F„  =ftL  where  0  is  the  number  of 
0K  molecules  incorporated  into  a  unit  volume  of  Si0B  and  L  is  the 
oxidation  rate.  The  diffusive  flux  of  0«  in  SiOs  has  not  been  measured 
directly.  However,  F0  ie  related  to  the  parabolic  rate  constant 
defined  in  the  Deal-Grove  model  <14>.  Thus,  an  approximate  value  was 
calculated  using  previously  published  kp  values  at  lOOO^C  and 
extrapolated  values  at  700°C  (9) .  The  kp  value  used  for  llOO^C  was 
determined  from  the  data  collected  in  this  experiment.  The  calculated 
ratio  Fr>  /  F«  in  several  thickness  regimes  is  shown  in  Table  3  for  the 
(100)  orientation  at  TOO^C  and  1000=C  and  for  the  (111)  orientation  at 
HOO^C.  Similar  results  were  obtained  on  the  (110)  orientation.  In  all 
cases,  the  ratio  FD  /  F«  is  greater  than  one  which  means  that  the 
oxidation  ia  not  diffusion  limited  although  a  diffusion  limited  regime 
is  approached  at  HOO'^C  for  thicknesses  above  300  nm.  The  lack  of  a 
direct  measurement  of  the  diffusivity  and  the  validity  of  applying  a 
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Fickian  diffusive  mechanism  to  the  thin  films  obtained  at  700“C  are 
sources  of  error  in  the  calculation.  However,  based  on  the  available 
data  used  for  this  calculation,  it  is  concluded  that  the  oxidation 
rate  is  controlled  by  the  interface  reaction  rate  in  the  thickness 
regimes  covered  by  these  experiments. 

To  determine  whether  the  Si  surface  atom  density,  N«i,  was  an 
important  parameter  controlling  the  oxidation  rates,  the  order  of  the 
rates  was  compared  with  NSi  as  a  function  of  orientation.  If  the 
interface  reaction  rate  is  dependent  on  the  Si  surface  atom  density, 
Nei,  (25)  then  the  experimental  rates  on  the  five  orientations  should 
parallel  the  order  of  Nai  with  orientation.  Calculated  values  for  NBi 
as  a  function  of  orientation  were  published  by  Ligenza  (1)  and  are 
shown  in  the  second  column  of  Table  4  along  with  the  value  for  Nai  on 
the  (511)  plane  which  was  calculated  using  the  same  method.  A  small 
correction  to  the  previously  published  Nei  value  for  the  (311) 
plane  (1)  has  been  made,  but  does  not  change  the  order  of  NSi  with 
orientation.  The  method  used  to  calculate  these  numbers  is  described 
in  Appendix  1.  It  is  seen  from  column  2  of  Table  4  that  the  order  of 
the  Nsi  values  calculated  using  this  method  is: 

(511)  >  (311)  >  (110)  >  (111)  >  (100) 

which  differs  from  the  experimental  order  of  the  oxidation  rates  (see 
Table  1).  The  apparent  discrepancy  between  the  calculated  and 
experimental  results  can  be  explained  using  the  results  of  LEED  and 
RHEED  studies  on  the  (311)  and  (511)  planes  (35,36).  These  planes  are 


vicinal  aurfacas  to  tha  (100)  plana  and  are  not  atructurally  atable. 
The  Oil)  and  (511)  aurfacea  reconatruct  to  fora  a  atapped  surface 
having  (100)  oriented  terracea  and  (ill)  oriented  atep  edgea.  The 
areal  ratio  of  (111)  to  (100)  oriented  regiona  on  the  reconatructed 
surface  ia  determined  by  the  initial  interplanar  angle  between  the 
vicinal  surface  and  the  (100)  plane.  Since  values  for  N«i  on  the 
reconatructed  (311)  and  (511)  surfaces  have  not  been  reported,  a 
calculation  was  made  based  on  the  available  information  (see 
Appendix  2).  The  results  are  shown  in  Column  3  of  Table  4.  These 
calculated  values  represent  effective  surface  atom  densities  since 
there  are  actually  regions  of  two  different  orientations  on  the  same 
surface.  Using  the  corrected  values  for  Nw,  on  the  reconstructed  (311) 
and  (511)  surfaces,  qualitative  agreement  is  obtained  between  the 
order  of  the  initial  rates  and  the  Si  surface  atom  density.  The 
initial  rates  on  the  five  orientations  increase  in  the  same  order  as 
the  actual  Si  surface  atom  density. 

Table  5  presents  a  quantitative  analysis  of  the  relationship 
between  the  surface  atom  density  and  the  relative  rates  in  several 
thickness  ranges  at  700^C  and  1000,=>C.  All  ratios  are  with  respect  to 
the  (100)  surface.  It  is  clear  from  Table  5  that  although  the  initial 
oxidation  rates  at  700,r,C  scale  qualitatively  with  N«»,  they  do  not 
scale  quantitatively.  The  rate  ratios  are  substantially  greater  than 
predicted  based  on  Nm  which  indicates  that  surface  atom  density  is 
not  the  only  parameter  affecting  the  initial  oxidation  rate.  The  lack 
of  a  quantitative  scaling  may  reflect  differences  in  reaction  rates 
caused  by  the  orientation  of  bonds  on  the  surface,  stress  effects  or 


other  properties  ralatad  to  the  raactivity  of  tha  aurfaca  atoaa.  Tha 
affact  of  stress,  if  significant,  ahould  ba  graataat  at  tha  lowar 
oxidation  taaperaturaa  whara  tha  atraaa  ia  largaat.  Howavar,  baaed  on 
the  raaulta  at  700°C  and  1000°C,  it  ia  alao  clear  that  atraaa  does  not 
affect  the  order  of  the  ratea  before  the  croaaover.  The  (110),  (311) 
and  (100)  orientationa  have  similarly  high  intrinaic  atreaa  values 
while  the  (111)  orientation  has  the  loweat  measured  stress  (11,24). 

The  order  of  the  initial  ratea  does  not  correlate  with  this  trend  in 
stress  and  is  not  changed  when  stress  levels  are  reduced  at  higher 
temperatures.  Thus,  while  stress  may  alter  the  absolute  rates,  the 
overall  order  of  the  initial  rates  remains  unaffected.  It  is  concluded 
that  the  surface  atom  density  is  the  most  important  parameter 
establishing  the  order  of  the  initial  rates. 


Extending  the  analysis  to  thicker  films  using  the  1000“C  data,  it 
is  seen  from  Table  5  that  beyond  SO  nm  there  is  relatively  good 
agreement  between  the  rate  ratios  and  the  ratio  of  Net  on  the  (110), 
(311),  (511)  and  (100)  surfaces.  The  rate  ratio  on  the  (111)  surface 

is  higher  than  predicted.  Thus,  except  for  the  (111)  orientation,  it 
appears  that  even  beyond  the  crossover  thickness  the  Si  surface  atom 
density  is  important  in  controlling  the  oxidation  rates  for 
thicknesses  up  to  200  nm.  The  anomolously  high  rate  on  the  (111)  is 
related  to  the  crossover  which  is  discussed  below. 


The  crossover  reflects  a  thickness  dependent  parameter  which 
either  increases  the  rate  on  the  (111)  surface  relative  to  the  (110), 
or  decreases  the  rate  on  the  (110)  surface  relative  to  the  (111).  The 
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Mechanism  causing  tha  crossover  is  not  well  understood.  To  date  the 


only  model  for  this  behavior  involved  a  stress  enhanced  interface 
reaction  on  the  (111)  orientation  (10),  but  recent  stress  measurements 
showed  that  the  (110)  orientation  has  a  higher  intrinsic  stress  than 
the  (111)  orientation  (11,24).  Thus,  the  crossover  cannot  be  caused  by 
a  streaa  enhanced  interface  reaction  on  the  (111)  surface.  However, 
the  results  from  the  llOO^C  oxidations  reported  here  indicate  that  the 
crossover  may  be  related  to  the  intrinsic  stress.  At  1100*»C,  where 
intrinsic  stress  is  reduced  to  its  lowest  value  (11,24,25),  a 
crossover  is  not  observed.  Since  the  crossover  observed  at  lower 
temperatures  is  not  due  to  a  stress  enhanced  reaction  on  the  (111),  an 
alternate  explanation  based  on  a  stress  related  decrease  in  rate  on 
the  (110)  surface  is  proposed. 

The  decrease  in  rate  on  the  (110)  aurface  ia  attributed  to  a 
decrease  in  the  diffusion  rate  through  the  SiOa  film  since  a 
compressive  oxide  stress  is  expected  to  lower  the  diffusion  rate 
(27,37).  It  was  determined  above  that  the  oxidation  rate  is  controlled 
by  the  interface  reaction  in  the  thickness  range  considered  here.  A 
simple  first  order  rate  law  for  this  process  is  given  by  F*  =  k,X, 
where  k.~  is  the  interface  reaction  rate  constant  and  is  the 
concentration  of  oxidant  at  the  interface  (14).  Although  it  was 
concluded  that  the  initial  order  of  the  rates  is  determined  by  the 
surface  atom  density  (through  k,-),  Ct  may  be  affected  by  stress  as  the 
oxide  thickness  increases.  If  Ci  is  reduced  due  to  a  decreased  rate  of 
diffusion,  then  a  corresponding  decrease  in  the  interface  reaction 
rate  is  anticipated.  This  is  equivalent  to  reducing  the  partial 
pressure  of  0*  in  the  ambient.  Since  the  (110)  has  the  highest 


intrinsic  strsss  (10,24),  the  diffusion  rats  and  thus  Ci  would  bs 


rsducsd  ths  aost  on  this  oriantation.  Tha  diffusion  procass,  and 
tharafora  C*,  on  tha  (111)  oriantation  would  be  least  affected  by 
stress.  If  Ci  on  the  (111)  oriantation  becomes  higher  than  on  the 
(110)  orientation,  then  this  difference  may  compensate  for  tha 
difference  in  NSi  and  result  in  the  crossover  in  rate  to  (111)  > 

(110).  This  model  assumes  that  the  stress  level  is  sufficient  to  alter 
the  diffusivity.  A  higher  diffuaivity  on  the  (111)  orientation 
relative  to  the  (110)  and  (100)  orientations,  as  required  by  this 
explanation,  is  supported  by  kp  values  derived  from  applying  the  Deal- 
Grove  model  to  oxidation  data  in  the  temperature  range  800-1000,:>C  (9)  . 
In  the  Deal-Grove  model  kp  values  are  directly  proportional  to  the 
diffusivity  of  0*  through  the  oxide.  kp  values  determined  at 
temperatures  below  950°C  (9)  are  the  highest  on  the  (111)  orientation 
and  thus  indicate  that  at  lower  temperatures  the  diffusivity  of 
through  the  oxide  on  the  (111)  is  greater  than  on  the  (110) 
orientation.  Although  this  model  for  the  crossover  will  require 
further  experimental  justification,  it  is  consistent  with  the 
currently  available  data. 

SUMMARY 


It  is  concluded  that  the  Si  surface  atom  density  is  the  most 
important  surface  property  establishing  the  order  of  the  initial 
oxidation  rates.  A  quantitative  scaling  of  oxidation  rate  with  surface 
atom  density  is  not  observed  in  the  thickness  range  0-30  nm.  Thus, 
other  orientation  dependent  parameters  alter  the  absolute  rates  but  do 


not  change  the  order  of  the  rates  in  the  thickness  regime  before  the 
crossover.  For  film  thicknesses  beyond  the  crossover  at  1000*»C  the 
ratea  scaled  in  approximately  the  same  ratios  as  the  surface  atom 
density  except  on  the  (111)  orientation.  An  explanation  for  the 
crossover  in  rate  between  the  (110)  and  (111)  orientations  is  still 
not  clear.  However,  at  1100°C  a  crossover  is  not  observed  and  the 
initial  order  of  the  rates,  (110)  >  (111)  is  maintained  at  thicknesses 
up  to  700  nm.  A  model  for  the  crossover  based  on  a  stress  related 
decrease  in  the  diffueivity  on  the  (110)  plane  at  lower  temperatures 
was  proposed. 


Appendix  1 


Calculating  the  Si  surface  atom  density  using  the  method 
described  by  Ligenza  (1)  involves  determining  the  number  of  Si  atoms 
located  on  a  specific  plane  in  the  Si  unit  cell,  and  dividing  by  the 
area  of  that  plane.  The  area  of  the  (100)  and  (110)  planes  is  easily 
determined.  The  calculation  for  the  area  of  the  (311)  plane  is 
described  below  since  the  value  determined  here  differs  from  Ligenza's 
value  (1).  The  area  of  the  (511)  and  (111)  planes  directly  parallels 
the  calculation  on  the  (311)  plane.  Referring  to  Figure  2  where  a  is 
the  lattice  parameter  of  the  Si  unit  cell,  the  area  of  the  (311)  plane 
is  calculated  as  follows: 

AB  =  BC  =  »f(  a®  ♦  (a/3)“  >  =  ^(10>*a/3 
AC  =  »T(2)*a 

BD  =  ✓(  (AB)*  -  (AD)®  )  *  tf(  10/9*a  -  l/2»a) 

So, 

BD  =  l/3*J-(ll/2)  *a 

Area  =  0 . 5* ( AC ) * ( BD )  *  1 /6*V ( 1 1 ) *a® 

The  areas  of  all  five  orientations  are  given  in  Table  6. 

The  next  step  is  to  determine  the  fractional  number  of  Si  atoms 
located  on  each  plane  within  one  unit  cell.  The  calculated  value 
includes  only  the  circular  projection  of  the  Si  atom  (assumed 
spherical)  on  the  plane  area.  For  the  (100),  (110)  and  (111)  planes 

the  calculation  is  direct,  but  for  the  (311)  and  (511)  orientations 
there  is  some  ambiguity  in  determining  the  number  of  Si  atoms  on  the 
plane.  On  both  planes  the  atom  located  at  (0,X,%),  D  in  Figure  2, 


contributes  O.S  to  the  totel  nuaber  of  Si  atoms.  For  the  (311)  plane, 

the  atoms  located  at  (0,1,0)  and  (0,0,1),  A  and  C  in  Figure  2, 

contribute  the  fraction  47.9/360  to  the  number  of  atoms  on  the  plane 
since  the  angles  ACB  and  CAB  are  47.9°.  On  the  (511)  plane  these  atoms 
contribute  the  fraction  46.1/360.  The  ambiguity  in  counting  the  number 
of  Si  atoms  on  these  surfaces  arises  from  the  treatment  of  the  atom 
located  at  (.%,%,%'>.  This  atom  does  not  lie  directly  in  the  plane  for 
these  two  orientations  but  is  very  close.  The  perpendicular  distance 
from  the  atom  at  (%,%,%)  to  the  (311)  plane  is  0.04  nm  and  to  the 
(511)  plane  is  0.06  nm.  In  calculating  the  surface  atom  densities  on 
these  two  plane  it  was  assumed  that  this  atom  contributed  1  to  the 
total  number  of  Si  atoms  even  though  it  is  not  directly  on  the  plane. 
Thus,  the  number  of  atoms  on  the  (311)  plane  is  1.77  and  on  the  (511) 

plane  is  1.76.  As  discussed  previously  in  the  paper,  the  results  of 

LEED  and  RHEED  studies  on  these  surfaces  (35,36)  show  that  there 
really  is  no  ambiguity  concerning  the  contribution  of  this  atom.  The 
structure  of  the  as  cleaved  surfaces  on  the  (311)  and  (511)  planes  are 
not  stable  and  reconstruct  to  form  a  surface  consisting  of  (100)  and 
(111)  regions.  A  summary  of  the  calculated  surface  atom  densities  not 
considering  any  changes  to  the  surface  structure  upon  cleaving  is 
presented  in  Table  6.  The  corrected  values  for  on  these  two  planes 

taking  into  account  the  reconstructed  surface  are  calculated  in 


Appendix  2 


In  order  to  calculate  the  effective  surface  atom  density  on  the 
(311)  and  (511)  planes  it  is  necessary  to  calculate  the  areal  ratio  of 
(100)  to  (111)  oriented  regions  on  the  reconstructed  surfaces.  This 
calculation  is  based  on  descriptions  of  the  (100)  vicinal  surfaces 
provided  by  LEED  and  RHEED  studies  (35,36).  Figures  3  and  4  depict  the 
essential  features  of  the  calculation  for  the  (311)  plane.  Consider 
the  area  ABCD  in  Figure  3  which  has  the  (311)  orientation  and  extends 
over  two  unit  volumes,  w3,  of  the  Si  crystal  where  w  is  an  integral 
number  of  unit  cell  lengths.  The  edges  of  the  unit  volume  are  oriented 
along  <100>  directions.  According  to  the  model  for  the  vicinal  (100) 
surfaces  (35),  the  original  (311)  plane  reconstructs  to  form  a  more 
stable  surface  consisting  of  (100)  terraces  and  (111)  steps.  This 
surface  structure  is  represented  by  the  dashed  lines  in  Figure  3  where 
a«ioo>  represents  the  area  of  each  terrace  and  a,,,,,  is  the  area  of  a 
step.  The  ratio  of  the  area  of  (100)  to  (111)  oriented  regions  in  ABCD 
is  given  by 

A  <  i  o*c>  >  /  A  c  i  i  i  >  =  Z  a  <  x  o  o  >  /  I  a  <  i  i  i  >  (1). 

The  areas  a<ioo>  and  aclii>  are  parallelograms  having  the  same  height, 
h,  but  different  widths  w<]Lo<>>  and  w<,,,>.  The  area  of  each  terrace  is 
then 

a<ioo>  =w<ioo>  xh  (2a) 


and  the  area  of  a  step  is 

a , i t i ,  =  w  < i i i >  xh  (2b). 
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A  <  i oo >  /  Amt)  *  Ituooi  /  ZSdU)  *  Iwtoo  /  Zwiii  (3) 

and  the  problem  reduces  to  the  2  dimensional  situation  depicted  in 
Figure  4  which  shows  the  region  OCD  from  Figure  3.  The  angle  8  is  the 
angle  between  the  <110>  direction  (i.e.  the  intercept  of  the  (111) 
plane  with  the  <001)  plane)  and  <010>  direction  which  is  45°.  It  is 
necessary  to  calculate  the  ratio  £w<ioo>  /  Ew(lll)  under  the 
reatrictions  : 

Ih.t.p  =  w/3  <  4  > 

and 

Iw < i oo >  *  Ew ' < i i i )  =  w  (5) 

where  h.«.p  is  the  step  height,  w/3  is  the  distance  0C  where  the  (311) 
plane  intersects  an  edge  of  the  unit  volume  and  w'(111>  is  the 
projection  of  the  (111)  step  edge  along  the  <100>  direction.  The 
number  of  steps  on  this  surface,  N_ ,  is  given  by 
N.  =  (w/3)  /  h.^.p  (6)  . 

Other  useful  expressions  from  figure  4  are  : 
w  <iii)  =  h .  *  „ /  tan(45°)  =  h.*  (7) 

and 

w  <  1 1  i  >  =  h«t /  sin(45,=  )  =  J’2»h»«.p  (8). 

Using  the  above  expressions,  £w<10o>  and  Ew , , , , >  can  be  derived  as 
follows : 

Rearranging  equation  (5), 

Ew (  \ oo  >  ~  W  —  Iw  i  x  x  x  >  (9). 

Using  equations  (6)  and  (7), 

Ew'  uui  =  N.  x  w'c*i»>  =  w/3  so, 

£w,ioo>  =  w  -  w/3  =  2  •  w  /  3  (10). 

Using  equations  (6)  and  (8), 


(2  C) 
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Thus,  the  areal  ratio  of  (100)  to  (111)  oriented  regions  on  the 
reconstructed  (311)  surface  is 

Zw  <  1  OO  >  /  Zw  <  i  i  i  >  =  (  2»w/3 )  /  (*/2»w/3)  =  1.414. 

For  the  (511)  plane  the  calculation  is  identical  to  that  on  the  (311) 
plane  except  the  number  of  steps,  N«  is  (w/5)/h.*.p,  instead  of 
(w/3) /h«t-r>  •  The  areal  ratio  (100):  (111)  on  the  reconstructed  (511) 
surface  is  2.83.  An  effective  Si  surface  atom  density,  Nai>  can  be 
calculated  for  these  surfaces  using  the  areal  ratios  and  the  Net  on 
the  (100)  and  (111)  planes: 

Net  -  A  (  i  iji  *  >  (  A  t  ot  *  i  <  i  utj  >  +  A  t  i  i  i  j  /  Atot  *  Net  till) 

where  At  .tt  ~  A  1 1  oo >  *  A  <  t  t  i  >  . 

For  the  (311)  orientation  this  is  calculated  to  be  7.21  x  lO1* 
atoms/cm^  and  for  the  (511)  orientation  7.05  x  101"*  ato»s/c»e. 
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Tabla  1 


OXIDATION  RATES  AT  700“C  in  (na/ain>  AS  A 
FUNCTION  OF  OXIDE  THICKNESS  FOR  FIVE 
Si  ORIENTATIONS. 


OXIDATION 

RATE  AT  THE 

FOLLOWING 

THICKNESSES: 

SILICON 

ORIENTATION 

5  na 

10  na 

20  na 

30  na 

(110) 

0.196 

0.129 

0.094 

0.074 

(111) 

0.130 

0.102 

0.088 

0.077 

(311 ) 

0.125 

0.089 

0.066 

(511) 

0.103 

0.074 

0.058 

(100) 

0.075 

0.057 

0.044 

TABLE  2. 

PARTIAL  THICKNESS  VERSUS 
OXIDATION  AT  1100~C. 

TIME  DATA 

SET  FOR 

OXIDE 

THICKNESS 

(na) 

OXIDATION 

TIME  (HRS)  (110)  (111) 


Tabla  3 


RATIO  OF  DIFFUSIVE  FLUX  TO  EXPERIMENTAL  FLUX  IN 
SEVERAL  THICKNESS  REGIMES. 


X0.  <na) 


Fd  /  F.kp 


(100)  Si  at  700°C 

5 

7.44 

10 

4.90 

20 

3.18 

<100)  Si  at  1000°C 


30 

50 

200 


9.29 

6.81 

4.70 


(111)  Si  at  1100~C 

100 

1.49 

300 

1 . 16 

500 

1 . 10 

700 

1.07 

TABLE  4.  Si 

SURFACE  ATOM  DENSITV 

VERSUS  ORIENTATION 

N*i  i  BY  METHOD 
OF  LIGENZA  <1) 

< APPENDIX  1) 

CORRECTED  N„* 
BASED  ON  LEED 
MODEL  <35,36) 

< APPENDIX  2) 

SILICON 

ORIENTATION 

Si  atons/cn* 
<10**) 

Si  atoms/cm® 
<101"*) 

<110) 

9.59 

<111) 

7.83 

<311) 

10.86 

7.21 

<511) 

11.48 

7.05 

<100) 

6.78 

TABLE  5.  A  QUANTITATIVE  COMPARISON  OF  OXIDATION  RATE  AT  700-C 
AND  1000~C  WITH  Si  SURFACE  ATOM  DENSITY 


ORIEN¬ 

TATION 

RATE  RATIOS  RELATIVE  TO  THE  <100) 

AT  THE  FOLLOWING  OXIDE  THICKNESSES: 

Nai  RELATIVE 
TO  THE  <100) 

700~C 

1000*=*C 

10  nm 

20  nm 

50  nm 

300  nm 

<110) 

2.26 

2.14 

1.51 

1.22 

1.41 

<  111 ) 

1.79 

2.00 

1.60 

1.45 

1.15 

<311 ) 

1.56 

1.50 

1.20 

1 . 14 

1.06 

<511) 

1 .30 

1.31 

1.08 

1.07 

1.04 

<  100) 

1 .00 

1.00 

1.00 

1.00 

1.00 

Table 

6.  CALCULATED 

SILICON  SURFACE  ATOM 

DENSITY 

ORIEN¬ 

TATION 

PLANE  AREA  in 

UNIT  CELL 

ATOMS  ON  PLANE 

IN  UNIT  CELL 

Si  atons/cm* 
(10**-) 

<110) 

V2a- 

4.0 

9.59 

<111) 

1  /2V3ac" 

2.0 

7.83 

<311  ) 

l/GVlla-* 

1.77 

10.86 

<511 ) 

3/10V3a- 

1.76 

11.48 

<100) 

2.0 

6.78 

a  =  5.431 

x  10"**  cm 

LIST  OF  FIGURES 


1)  Figura  1 


Oxide  thickness  (nm)  vtraui  oxidation  tiaa  (hours) 
for  700°C  silicon  oxidation  in  dry  oxygan.  Data  is 
shown  for  fiva  silicon  oriantations . 


2)  Figura  2. 


Tha  (311)  plana  within  the  silicon  unit  call  with 
lattice  parameter  a. 


3)  Figura  3. 


Tha  reconstructed  (311)  surface.  Area  ABCD  has  tha 
(311)  orientation  within  two  unit  volumes,  w3,  of  the 
Si  crystal  where  w  is  an  integral  number  of  unit  cells. 
The  dashed  lines  represent  the  reconstructed  surface 
containing  (100)  terraces  and  (111)  steps  (35,36) . 


4)  Figure  4. 


A  two  dimensional  representation  of  the  reconstructed 
(311)  surface  shown  in  Figure  3.  The  height  of  each 
step  is  h«o„  and  the  angle  e  is  45~. 
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